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Abstract 

In this paper I consider sharp spectral asymptotics for multidimensional magnetic 
Schrodinger operator with irregular coefficients with respect to two parameters - semi- 
classical parameter h and coupling parameter ^. There are few principally different 
cases, depending on dimension, rank of magnetic intensity matrix, relation between 
h and /i and some extra assumptions. 

Introduction 
0.1 Preface 

In this paper I consider multidimensional Schrodinger operator 
(0.1) A = Ao + V{x), Ao=Y, P39'\x)Pk, 

j,k<d 

Pj = hDj - /iV,(x), h e (0, 1], /i > 1. 

It is characterized by magnetic field intensity matrix (Fj^) with Fj^ = (<9x Vfc— S^^^. V^), which 
is skew-symmetric d x d-matrix, and (F^) = {g-'''){Fkm) which is unitarily equivalent to 
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skew-symmetric matrix [g^^)^ [Fjk)[g^^)^ . Then all eigenvalues of (F^) (with multiplicities) 
are ii/m, (/m > 0, m = 1, . . . , r) and of multiplicity q = d — 2r where 2r = rank(F^). 

I formulate results only in the case of g^^ = 6jk (thus covering the case g^^ = const as 
well) and Fj^ = const (then one can select linear vector potential (Vi(x), . . . , Vd{x))). The 
results in the general case differ even in their statements. So we consider operator 

(0.2) J2 + '^(^)' = - l^^j(^) 

with linear functions Vj{x) and constant magnetic intensity matrix Fj^. 

The main goal of this paper is to present the local spectral asymptotics, i.e. asymptotics 

of 

(0.3) / e{x,x,T)il){x) dx 



as h ^ +0, fi — > +OC where e{x,y,T) = eh,i_i{x,y,T) is the Schwartz kernel of the spectral 
projector of A {A is assumed to be a self-adjoint operator) and ip is a smooth cut-off function 
supported in the ball B{0, |) (all the conjectures save self-adjointness are made for B{0, 1)). 
Combined with partition-rescaling arguments such asymptotics imply many asymptotics of 
eigenvalue counting function given by formula (0.3) with ip = 1. 

Under above assumptions in the appropriate coordinates after gauge transform 



(0.4) V,ix) 



— fj-j-Xj-r as J — r -|- 1, . . . , 2r 

otherwise 



where fj > and ±«/j are eigenvalues of F^, rankF = 2r. Let q = d — 2r. 

In the special case V = const everything becomes explicit: after /i- Fourier transform 
with respect to (x^+i, . . . , Xd) and change of variables Xj i— > Xj + fJ'~^ fj~^^2j+i {j = ^, ■ ■ ■ ,r) 
operator A is transformed into 

(0.5) Yl (^'^1 + + E + ^- 

l<i<r l<k<g 

Using Hermite function decomposition one can calculate easily e{x, y, r) and prove that 
(0.6) e{x,x,r)=S^:^{r)'^ 



UJq 



(27r)-^ E - + ^)fjf^h -vf X ifh~'^^h ■■■frV9 
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where Uq is a volume of the unit ball in M'^, g = det{g^^)~^. 

In particular, for q = such operator has pure point spectrum of infinite multiplicity, 
consisting of Landau levels = ^j(2aj + l)fjfih + V with a G Z"^*^' while for g > 1 the 
spectrum is absolutely continuous and Landau levels are merely bottoms of its channels 

Also, as fih ^ 1 one needs to include — Ylj fjf^h or even larger negative number in V 
to avoid being in the classically forbidden zone. 

Our goal is to study operator A with variable potential V{x); however the deep distinc- 
tions between cases rankF^fc = d and rankFj^ < d and between cases ///i -C 1, /x/i ^ 1 
(and intermediate case /i/i ~ 1) are preserved as well as the importance of Landau levels. 
Also, gives a good approximation to e(x, x, r) even if in some cases we need some 

corrections to derive sharper remainder estimate. 

0.2 Classical Dynamics 

Considering classical trajectories on some finite energy level one can prove easily that for 
V = const 

(i) As rankF = d = 2 particles move along circles of the radii x /i"^. As rankF = d = 
2r, r > 2 trajectories are more complicated (generic trajectories are periodic if fi, . . . , fr 
are commensurable or envelope tori otherwise) but they are confined to C/i~^-vicinities of 
their origins as well; 

(ii) As d = 3, rankF = 2 particles move along spirals. Similar description holds for 
rankF = 2r < d, r > 1: there is a cyclotronic movement as in (i) along [xi, . . . ,X2r) and 
also a free movement along x'" = {x2r+i , • • • , Xd) ■ 

So, we see a big difference between full- and not-full-rank cases and more subtle difference 
between d = 2,3 and d > 4. 
As V is variable 

(iii) As rankF = d = 2 particles move approximately along circles of the radii x /i^^ 
but the centers of these circles are drifting with the speed /i~^|/i|~^VV^"'" (thus the drift 
is orthogonal to the direction of the electric field). Similarly, as rankF = d = 2r the 
trajectories described in (i) are drifting with the velocity fi~^F"^W where F = (Fjk) is 
the intensity matrix; 

(iv) As d = 3, rankF = 2 there is a fast cyclotronic movement along (xi, X2), temperate 
movement described by 1-dimensional Hamiltonian -|- V{x', X3) + -Emagn along ^3 where 
-E-magn is an energy of the fast movement (which is constant in our assumptions) and a slow 
drift along all variables. Similar description holds for rankF = 2r < d, r > 1. 

So, if rankF = d we can follow trajectories until time Ti = e/i. Then assuming that 
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trajectories are non-periodic which happens if 

(0.7) \V\ > eo, 

(0.8) |V\^| > eo 



we hope to get remainder estimate which is Tf h = n~^h times better than the principal 
part; since the latter is x h~'^(^l+fj,^h^^ we expect to derive asymptotics with O (/i~^/i^~'^(l + 
fi^h'^)) remainder estimate under assumptions (0.7)-(0.8). Assuming sufficient smoothness 
(see below) it is a correct guess. 

However without condition (0.8) trajectories could be periodic with period Tq = e/i~^ 
and we cannot expect remainder estimate better than 0(/i/i^~'^) as fih < 1. 

On the other hand, if rankF < d we can follow trajectories until time Ti = e and under 
assumptions (0.7)-(0.8) these trajectories are non-periodic. So, we expect 0(^h^~'^{l+fi'^h^)^ 
remainder estimate then. 

However, even without condition (0.8) most of the trajectories are non-periodic in a short 
run due to the free movement. Actually, in the classical settings periodic trajectories must 
have ^2r+i = ■ ■ ■ = = and thus form a set of measure but we are in semiclassics and 
talking about non-periodicity one should be able to observe the shift after "quasi-period" 
X Anyway, we guess that the lack of condition (0.8) would not be that destructive 

(and even noticeable for not very large fi). Both guesses are correct under some smoothness 
assumption. 

As /i/i -C 1 one can rescale to the standard case (i.e. with /i = 1). Really, after scaling 
jjix ^ X and thus hD i— > jjihD, i— > 1, /i yu/i we get a standard case and scaling back we 
get the Weyl principal part 

(0.9) X, r) = uJd{T - Vyjh-''^ 

and the remainder estimate 0{fih^~'^). Much more useful however is a rescaling applied to 
the intermediate rather than the final results. 

0.3 Canonical form 

Starting again from appropriate coordinates and applying /i-Fourier transform with respect 
to x" = (xr+i, ■ ■ ■ , Xr) and change of variables xj i— > xj + fi~^f~^^2j+i with respect to 
x' = (xi, . . . , Xj.) we arrive instead of (0.5) to 

(0.10) (^'^1 + f^'fM) + E ihD2r^^?+ 

l<j<r l<k<q 

V{x' + fi-^hP-^D", x" - fi''hF~^D', x'") 
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where the third term is /i~^/i-pseudo-differential operator. 

Note that on bounded energy levels operators hD' and fix' are now bounded. Then one 
can apply Taylor decomposition to the last operator in (0.10) leading to an operator with 
the main part 

(0.11) A'= J2 (h'D^ + t^'fM) + E ihD2r+,f + V{f,-'hF-'D",x",x"'), 

l<j<r l<k<q 

junior terms 

(0.12) A' = Yl Wai3{fi-'hF~^D",x",x'")x{fi~'hf\x"'D'^ 

a,PeZ+ r ,l<\a\ + \(3\<l 

and the remainder estimate 0{^~^) (subject to the smoothness assumptions). 

Applying Hermite decomposition with respect to x' we see that as g = the main part 
becomes a family of /i~^/i-pseudododifferential operators with respect to x" 

(0.13) Aa = Vc,= {2aj + l)fjijh + Vi^-^hF-^D", x" , x'"), 

l<k<r 

while as g > 1 it becomes a family of g-dimensional Schrodinger operators with respect to 
x'" with potentials which are yU~^/2,-pseudododifferential operators: 

(0.14) Aa= Yl ihD2r+kf + Vo,, 

l<k<q 

Va= Y {2aj + l)fjfih + V{fi-^hF-^D",x",x"'). 

l<k<rq 

As g = the principal part of the spectral asymptotics for an individual operator is 
X fi'^h"^ with the remainder estimate 0{fi'^~^h^~^) under non- degeneracy condition 

(0.15) \Va\ + \\/x",("'Va\>e 

and 0(/i''/i~^) otherwise (functions could be very flat). Assume that this non- degeneracy 
condition (which is equivalent to (0.8) as /x/i <C 1 and replaces it otherwise) holds. Then 
there are x [fi~'^h~^ + 1) contributing operators in the family (the rest is excluded due to 
ellipticity arguments) thus leading to the principal part and the remainder estimate x /i"*^ 
and 0(/i~^/i^~^^) respectively as yu/i < 1 and x fi'^'h~''' and 0{fi'^'~^h~^'^^) as fih > 1. 

For individual Schrodinger operator the principal part is x fi'''h^~'^ with the best possible 
remainder estimate 0{fi^'h^^'^~'^) under condition (0.15). However this condition is not 
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crucial: this remainder estimate holds without such condition if either g > 3 or g = 2, 
V ^ and one can recover a weaker remainder estimate otherwise. Again there are 
X [fi~'^'h~^ + 1) contributing operators in the family thus leading to the principal part 
and the remainder estimate x h~'^ and 0{h^~'^^) respectively as /i/i < 1 and x fi'''h~^' and 
0{n^h~^~^^) as fih > 1 (in the best case). 

Surely one can take care of junior terms (0.12) and remove them if possible, thus reducing 
operator to Birkhoff normal form. It is easy in the smooth case as d = 2, 3; as ci = 2 the 
normal form is 

(0.16) /i^ amjkix2,f-i-'hD2)-{xl + i2-'h^Dl)'^fi-'"'-'^-''h'' 

m+j+k>l 

and for = 3 it is 

(0.17) fi"^ ^ ampjk{x2,Xs,iJ~^hD2)x 

m+p+j+k>l 

X {xl + fi-'h^Diy\hD3)'y~'^~'^~^^~'h'' 
but as r > 2 the resonances become one of the obstacles: namely we cannot remove terms 
(0.18) a^^p{x", x'", n^^hD") x JJ {xj + ijj-^hDjp {xj - i^i"^hDjY^ x (^'")^, 

with 

(0.19) ^(«^._/5,.)/^. = 

i 

(c? = 3 is a special case). So resonance means that 7j/j = with 7 G Z*"; I7I is an order 
of resonance. 

As g = 0, the second order resonance terms are {xj + ifi~^C,j){xk — ifJ'~^^k) with fj = f^, 
the third order resonance terms are {xj+ifi~^^j){xk—ifi~^C,k){xm—ifJ'~^^m) with fj = fk+fm 
and their conjugates, etc. As g > 1 there are additional (but less malicious in the end of 
the day) resonance terms. 

Another obstacle is the lack of a very large smoothness. I overcome both of these 
obstacles because most of the junior terms are non-essential small perturbations: as g > 2 
all junior terms are non-essential; as g = 1 all terms but of the form ^~^{xj + i^~^^j){xk — 
ifJ^~^^k){xm — ifJ^~^^m) are are also non-essential (but such terms do not appear for special 
operator (0.2)), while for g = there are more essential terms. 

Further, the canonical form reduction is possible for general operator (0.1) as well but 
there will be essential cubic perturbations and third-order resonances become important as 
g = l. 
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0.4 Main tools and classification 

We will consider u{x,y,t) the Schwartz kernel of operator U{t) = e*^ where A is one 
of two framing approximations for A; in the case of irregular coefficients A will include 
mollification as well. However even in the case of smooth coefficients it is often convenient 
to take A different from A. 

Microlocal analysis (propagation of singularities) of u{x,y,t) in its original form or 
transformed (as A is reduced to its canonical form) plays a crucial role. The former approach 
is used in the case of weak magnetic field and it is combined with the standard theory 
rescaled. In the latter approach is used in the case of intermediate and stronger magnetic 
field and it is combined with the successive approximation method to construct solutions. 

As g = and non- degeneracy condition holds we need to take mollification parameter 
e = Cfih\ logh\ to make a weak magnetic field approach working properly. On the other 
hand, if we reduce operator to its canonical form, we must take e = C{fi~^h\\ogfi\)2 as 
d = 2 and e = C max(/i~^, {fi~^h\ log/i|)2) a.s d> 4. Since we are interested in the smallest 
possible e we apply weak magnetic field approach for /i < C{h\ \ogh\)~^ as d = 2 and for 
/i < C{h\\ogh\)~2 a.s d > 2. Otherwise we refer to intermediate magnetic field case as 
fi < C{h\ \ogh\)~^, then strong as /i < eh~^, superstrong as /i < Ch~^ and ultrastrong as 
/i > Ch'\ 

This classification holds as g > 1 but with a twist. If magnetic field is weak we study 
only "original" propagator u{x,y,t) picking up e = C p'^h\logh\ with p = as p > 
Pi = Cmax(pt~^, {fih\ \ogh\)2^ and e = Cp^^h\ \ogh\ otherwise. Inner zone ^ Pi} is 

treated by referring to the standard results rescaled. 

If magnetic field is intermediate we apply a weak magnetic field approach in the outer 
zone {|^'"| > Pi}- In the inner zone we apply reduction to the canonical form. Threshold 
between weak and intermediate magnetic field cases depends on q and the presence of 
non-degeneracy condition. 

In the strong (and stronger) magnetic field cases we use canonical form reduction ev- 
erywhere and e = Ch\ log/i|. 

Actually the above choice of e is needed to provide a reduction; the further analysis as 
q > 1 requires to increase e and the actual choice of e varies (as q = 1,2) depending on 
other assumptions and those indicated are kind of milestones rather than the actual values. 

1 Weak Magnetic Field 

We assume that magnetic field is relatively weak and we prove our results by reducing to 
the standard case. 
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1.1 Standard reference result rescaled 

As I mentioned the intermediate resuh in the standard case will be valuable for us more 
than the final one; namely we have by rescaling 

(1.1) \F,^,-^,{Mt)r.u)- Jn{iT-\)^)d,S'^'^ix,X)\<Ch'~'' Vr : |r| <e 

as /i < eh^^l log/i|~^, Ch\ logh\ < T < Tq = e/i^^; here and below x is smooth and even 
satisfying admissibility condition of [Brivr] function, supported in [—1, 1] and equal 1 in 
I]' X is its Fourier transform, V^v = v{x,x,t). 
We use also notations Tv = j T^ydx and (prit) = (f){t/T). Also we need to assume con- 
dition (0.7). Estimate (1.1) holds as coefficients of operator are either smooth or mollified 
with e > Ch\ logh\. 

1.2 Full-rank case 

We want to increase T in (1.1). To do this we analyze propagation of singularities. Let us 
introduce slow variables 

(1.2) Xj = Xj-^ /3jkPk{x, 0, Pk{x, = l^~^^k - Vk 

with {l3jk) inverse matrix to (Fj^). 
Then the Poisson brackets satisfy 

(1.3) {pk,Xj} = o yj,k, 

(1.4) {X„ Xk} = fx-'Pjk, {Pj,Pk} = -f^-'hF,k Vj, k. 

Then in the classical dynamics for time T the shift of 0(x, ^) = Ylj'^j-^j is equal 
to l^~^T f, PjkijdkV + o(l)) as T < Ti = e/i and under condition (0.8) it will be of 
magnitude fi~^T for an appropriate vector i, \i\ = 1. 

To make it observable from the point of view of microlocal analysis one needs to satisfy 
the logarithmic uncertainty principle (see [BrIvr]) which in this case is 

(1.5) fi~^T X e > Cfi~^h\ log h\ 

because Xi, . . . ,Xa are linked by (1.5) and semiclassical parameter is therefore fi~^h. We 
can avoid condition 

(1.6) exe>Cfx-^h\\ogh\ 
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because we do not use any reduction in this case; condition (1.6) be very unpleasant for 
not very large /i. 

We want to satisfy (1.6) with T = Tq and therefore we pick up the smallest e to satisfy 
this condition: 

(1.7) e = Cfih\\ogh\. 
Then we prove that 

(1.8) \Ft_.h^i,XT{t)r{u4j)\ < Ch' Vr : |r| < e VT G [fo,fi] 

where ip = ip{x) is a cut-off function described above, x is smooth (and [Brlvr] type) 
function supported in [— 1, — |] U [|, 1]. 

Then (1.1) and (1.8) imply that estimate 

(1.9) |F,^,-vXT(t)r(w^) - J J n{{r - \)L)d^£^'^{x, A) ij{x)dx\ < Ch"-' 

holds with any T G [To,Ti] and Vr : |r| < e. 

We must here and below integrate with respect to ip{x) dx because Xj contain and 
we need a real trace F and not just a restriction to the diagonal T^. 

Then Tauberian arguments imply spectral asymptotics with the remainder estimate 
0{T^^h^^'^) = 0{n~^h^~'^) while an approximation error is 0(e'| loge|~'^/i~'^) provided 

(1.10) VeC^'" {l,a)y (1,1) 

where C'''^ is a class of functions with k-th derivatives continuous with continuity modulus 
fi-k^ ^Qg^|-(T functions, k = [/J unless / G Z, a < when k = I — 1. 
Thus we arrive to our first statement: 

Theorem 1.1. Let d = 2r, V e C''^ with {I, a) h (1,2) and conditions (0.7), (0.8) hold. 
Then for 

(1.11) fi<h^~^ 

with an arbitrarily small exponent 6 > the following estimate holds 

(1.12) 7^ = I J (e(x, X, 0) - ^^W(x, 0))^{x) dx\ < 

Cfi-^h^-"^ + C{fih\ log/i|)'| log/i|-"/i-'^. 

In particular, as {l,cr) = (3,1) sharp remainder estimate TZ < Cfi~^h^~'^ holds for 
fi < h~2\logh\~2 and as = (2,1) this sharp remainder estimate holds for fi < 

h~3\ loe:/i|~3 . 
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1.3 Non- full-rank case 

As g > 1 we can use a free movement along x'" to extend (1.9) from T = Tq to larger T. Let 
us introduce p-partition with respect to C,'" with p = and let us consider a partition 

element with p > Cfi~^. Then for time T < Ti = ep the shift with respect to x'" will be of 
magnitude ep and in order to be observable it must satisfy logarithmic uncertainty principle 

(1.13) p X pT > Ch\\ogh\ ^ T>Cp~'^h\\ogh\. 
In order to be able to plug T = Tq = e/i~^ we must take 

(1.14) P>Pi = Cmax(/i-\(/i/i|log/i|)5). 

Then both (1.8) and (1.9) hold for uip replaced by Q{hDx){uip) with /i-pseudodifferential 
oprator Q with symbol supported in the indicated partition element. 

Then due to Tauberian arguments the contribution of this element to the remainder 
estimate is Ch^~'^T^^ p'^ and therefore the total contribution of the outer zone {|^'"| > pi} 
to the remainder estimate does not exceed Ch}~'^ f T^^{p)p'^~^ dp which is 0{h^~'^) as g > 2 
and 0{h^~'^\og fi) as g = 1. To get rid off the logarithmic factor we prove that for (/, a) >z 
(1, 2) one can take Ti = ep| logpp in the propagation. 

Also in this zone e must satisfy logarithmic uncertainty principle e x p > Ch\ logh\ and 
taking e = Cp~^h\ \ogh\ there we get mollification error Ch~'^ j \oge\^" p'^^^ dp which is 
0{h^-'^) as either g > 2, (/, a) h (1, 1) or g = 1, (/, 0") ^ (1, 2). Note that now mollification 
parameter e depends on (x, ^). 

On the other hand, contribution of the zone {|^'"| < Pi} to the remainder estimate is 
CT^^plh^-'^ = 0{h^~'^ + {ph)i+^\ \ogh\ih-'^). 

Finally, in this zone we pick e = Cp^^h\ logh\ and its contribution to the mollification 
error will be less than what we already got. 

So we arrive to estimate 

(1.15) 7^ < Ch^-"^ + C(p/i)2+i| log/i| U"-^ 
which actually could be improved to 

Theorem 1.2. Let condition {0.7} be fulfilled. Let either q = 1, V e C^-^ or q > 2, 
V e C^'\ Then 

(1.16) n<Ch^~^ + C{fih)i+^h-'^. 

In particular, for p < p^^^ = h 9+2 sharp remainder estimate TZ < Ch holds. As 
g > 2 this sharp remainder estimate holds for ji < . 
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Note, that at this stage extra smoothness is not very usefuL Furthermore we can assume 
in what follows that /i > /i^^-j and thus pi = C{fih\ log/i|)2. 

To improve estimates (1.15), (1.16) one needs to use better arguments in the inner zone 
^ = — Pi}- Under non- degeneracy condition (0.8) we can apply the same arguments 
as in the proof of theorem 1.1 as long as e = Cfih\ logh\ in Z. Then contribution of this 
zone to the remainder estimate becomes 0{h^~'^) while its contribution to an approximation 
error becomes 0{p\e^ \ \oge\~''h~'^) and we arrive to 

Theorem 1.3. Let condition (0.7), {0.8) he fulfilled. Let q>l, V & C^'"^ and either q = I, 
(/, a) y (1, 2)orq> 2, (/, a) >z (1, 1) . Then 

(1.17) 7^ < Ch^-'^ + C{fih)^+^\ log/i|'-"/i-'^. 

In particular, sharp remainder estimate TZ < Ch^^'^ holds for /i < h^^\\ogh\~^ as 
q=l, (/,a) = (|,i). 



2 Intermediate and strong magnetic field 

Now magnetic field is stronger than before but still either below eh~^ \ \ogh\~^ or between 
this value and eh~^. There is certain difference in the analysis of these two cases and for 
general operator (0.1) some statements would slightly differ as well. 



2.1 Full- rank case 

In this case we pick 

C(/i-^/i|log/i|)5 r = l, 

Cmax(/i-\ (p~^/i| log/i|)5) r>2 



(2.1) 



and reduce operator to a canonical form in the smooth case or "a poor man's canonical 
form" otherwise. Cases r = 1 and r > 2 differ because of the reduction: we need to solve 
equation 

(2.2) {ao, S} = V-W, ao= J2 /^ (s' + 

i<i<'- 

where W consists of unremovable terms; as r = 1 this equation is solved by integration 
along circles leaving W = W(^X2,C2,{x1 + ^1)2) while for r > 2 it is solved by Taylor 



2 INTERMEDIATE AND STRONG MAGNETIC FIELD 



12 



decomposition and we must assume that e > Cjj . On the other hand, we always need 
to assume (1.6) now because we consider x",^" as dual variable and we need to consider 
/i~^/i-pseudodifferential operators. 

In the best possible case we would get something similar the family of separate scalar 
yU~^/i-pseudodifferential operators but the same results hold in the general case as well: 

Theorem 2.1. Let d = 2r, V e C"'^ with {I, a) h (1,2) and conditions {0.7), {0.8) hold. 
Then 

(i) For fi < h~^\logh\~^ estimate 
(2.3) n <Cfi~^h^-'^ + C epilog hi'" h-'^ Vr : |r| < e 

holds with e = 

(a) For h~^\ \ogh\~^ < < eh~^ estimate {2.3) holds with e = C{n~^h\ log/i|)2 . 

In particular, as (/,cr) = (3, 1) sharp remainder estimate estimate TZ < Cfi~^h^~'^ holds 
for h~2 \ log/i|~2 < /i < eh^^ ; 

(Hi) As r = 1, estimate {2.3) holds with e = {fi~^h\ logh\)2 but for TZj instead ofTZ; 
here and below TZj is defined by formula {L12) but with S^"^ replaced by Sf^^ which is 
defined by {0.6) with V{x) replaced by W{x) where W{x) is an average ofV{y) along circle 
Cx = {y '■ \x — y\ = {fif)~^{T — V{x))2}; f is a scalar intensity of magnetic field. 

In particular, as {I, a) = (2, 1) sharp remainder estimate estimate IZj < C^~^h^~'^ holds 
for fi > h~3 \ log/;.|~3. 

This statement together with theorem 1.1 cover case /i < eh~^ completely. 
2.2 Non-full rank case. I 

The same classification and definition of e persist as g > 1; however we reduce operator 
to a canonical form only in inner zone Z = {\C,"'\ < Pi} (in the intermediate magnetic 
field case); in the outer zone we apply the weak magnetic field approach and pick up 
e = C\^"'\~^h\ logh\; furthermore after reduction is done, e is redefined (increased) in the 
inner zone as well. 

There are few different statements to prove; the first one is a generic one: 

Theorem 2.2. Let q > 1, V e C"'^. Then 

(i) As q > 3, (/, cr) = (1, 1) sharp remainder estimate IZ < Ch^~'^ holds for p, < Ch"^ ; 
(a) As q = 2, {I, a) = (1, 1) remainder estimate 



(2.4) 



7^ < ch^~^ + CiJ0-'^ 
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holds for jj, < Ch ; 

(in) As q = 1, {l,cr) >z (1,2) remainder estimate 

(2.5) 7^ < C^/i3-'^ + C/-5/ii-<i|log/i|-f 

holds for log/i|~^ < /i < Ch~^; 

(iv) As d = 3, (/, 0") = (1,2) remainder estimate 

(2.6) TZj < Cfihl"^ 

holds for < /U < Ch^^; here TZi is again defined by {1.12) with replaced by 
which is defined by {0.6) with V{x) replaced by W{x) where W{x) is an average of V{y) 
along circle Cx = {y : = ys, \x — y\ = {fif)^^{T — V{x))2} as magnetic field is directed 
along x^. 

In what follows we need to treat only cases g = 1,2. As non-degeneracy condition is 
fulfilled we get 

Theorem 2.3. Let q = 1,2, V e C^'" and conditions {0.7), {0.8) hold. Then 

(i) As q = 2, (/, cr) = (1, 1) sharp remainder estimate TZ < Ch^~'^ holds for fi < Ch~^; 
(a) As q = 1, (/,cr) y (1,2) remainder estimate 

(2.7) 7^ < C/i^^'^ + C/i^-'/i^-'^llog/il^"" 

holds for log/i|~^ < < Ch~^; 

(Hi) As d = 3, (/, 0") = (1, 2) sharp remainder estimate TZj < Ch^~'^ holds for fi < Ch~^ . 

So far all the the results of the article could be generalized to a general operator (0.1) 
with modification of non- degeneracy condition (0.8) and, as g = 0, 1, r > 2 with the special 
attention to the third order resonances (because in the general case they could lead to 
non-removable 0{n~^) terms in the canonical form). 

2.3 Non- full-rank case. II 

In this section we exploit more specific properties of operator (0.2), namely that fj have 
constant multiplicities. 

Theorem 2.4. Let q=l,2,Ve C^'" and condition {0.7) hold. Then 
(i) As q = 2, {I, a) >z (1, 1) remainder estimate 

(2.8) n<Ch^-'^ + CfihT7-2+^-'^\\ogh\-w-2 
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holds for /i < Ch~^; in particular sharp remainder estimate IZ < Ch^"'^ holds as {l^cr) = 
(2,0). 

(a) As q = 1, (1,2) ^ (IjCr) ^ (2,0) remainder estimate 

(2.9) 7^ < Ch^-'^ + Cfih^+^~'^\ log /ipiTs + log /i|-f 

/ioWs /or log/i|~5 < /i < Ch~^; 

(Hi) As d = 3, (1,2) ^ (/,cr) ^ (2,0) remainder estimate 

(2.10) 7^7 < C/i^-^ + C/i/i^+^-'^l log/iriT2 

/ioWs for I log < /i < Ch~^. 

Further, as r > 2 Diophantine properties of (/i, . . . , Z^) can play role. Assume that 

(2.11) n(n, r) #{a G Z+^ ^(2a, + 1)/, + V^^ < r} 

j 

satisfies estimate 

(2.12) \n{h,T) -n{h,T')\ < Ch~'' {\t - t'\ + u{h)) 

V/iG (0,1] Vr,r': |r| < C, |r'| <C 

with z/(?i) = o{h) (it holds with z/(/i) = /i for sure). Two following theorems improve 
theorems 2.2, 2.4 respectively: 

Theorem 2.5. Let r > 2, q = 1,2, V e C"'^ arwi conditions {0.7), {2.11) hold. Then 
(i) As q = 2, (/,cr) ^ (1, 1) remainder estimate 

(2.13) n<Ch^-'^ + Cu{iih)hi-'^ 

holds for /i < Ch^^ ; 

(a) As q = 1, {l,cr) y (1,2) remainder estimate 

(2.14) 7^ <C/i^-'^ + log/i|5--+ 

Cz/(/i/i) (^/i^-*^ + ^~h~'^\ log/i|~^) 

/ioWs /or /i < Ch~^ . 
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Theorem 2.6. Letr>2, q= 1,2, V G C'''" and conditions {0.7), {2.11) hold. Then 
(i) As q = 2, {I, a) >z (1, 1) remainder estimate 

(2.15) n < Ch^-'^ + Cv{^ih)h^~'^\ \ogh\-^ 

holds for n < Ch^^; 

(a) As q = 1, {l,cr) ^ (I72) remainder estimate 

(2.16) 7^ KCh^-"^ + log/i|5-'^+ 

Cv{^ih){hJ^~'^\ \ogh\~T^ + log/i|"?j 

/ioWs /or /i < Ch^^. 

Note that in estimates (2.13)-(2.16) right-hand expressions are sums of the right-hand 
expressions under non- degeneracy condition and of the right-hand expressions without mi- 
crohyperbohcity conditions, but latter are muhiphed by h'{fih){^h)~^. 

3 Superstrong and ultrastrong magnetic field 

In this case fi > eh~^ and the distance Landau levels increases; in the case of the ultrastrong 
magnetic field only one (may be multiple) level should be considered. 

3.1 Superstrong magnetic field 

In this case h~^e < fi < Ch~^ and the magnitude of the principal part of the asymptotics 
is still the same (/i"*^) as well as the remainder estimates in theorems 2.1-2.4. The only 
difference is that non- degeneracy condition (0.8) is replaced by 

(3.1) |r-l^-^(2aj- + l)/j/i/i| + |Vr| >eo Va G Z+^ 

j 

Theorem 3.1. (i) Statements of theorems 2.2, 2.4 remain true for h~^e < /i < Ch~^; 

(a) Statements of theorems 2.1, 2.3 remain true for h~^e < < Ch~^ with condition 
{0.8) replaced by {3.1) with r = and condition {0.7) skipped; 

(Hi) As q = under condition 

(3.2) \t-V -^{2aj + l)fji2h\>eo Va G 

j 

with r = estimate TZ < Ch"^ holds with arbitrarily large s (spectral gaps). 
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Remark 3.2. As 1 < fi < Ch ^ under condition 

(3.3) T-V-Y,fjf^h<-eo 

j 

estimate TZ < Ch'^ holds with £^^{x, r) = (r = 0) and with arbitrarily large s. 
3.2 Ultrastrong magnetic field 

In this case // > Ch~^ and in order not to be below the bottom of the spectrum one should 
modify condition V G C^'" . Assume instead that 

(3.4) V = - ^{2aj + l)fj^^h + W, a G W G C^'" {q = 0); 

j 

(3.5) y = - ^ /j/i/i + H/", H^gC'" (g>l). 

j 

Theorem 3.3. Let q = and conditions {3.4), {3.1) be fulfilled with r = 0. Then 
(i) Estimate 

(3.6) 7^ < c^'-^h-'+^- 

holds for fi > Ch^^ ; 

(a) Furthermore, under condition {3.2) estimate TZ < Cfi~^ holds with arbitrarily large 

s. 

Theorem 3.4. Let q>l and condition {3.5) be fulfilled. Then 

(i) Estimate 

(3.7) n < Ci/h-'+^ 

holds for fj, > Ch~^ , q > 3; 

(ii) Estimate 

(3.8) n < Ci/h-'+^{l + h-^+Th\ \ogh\-Th) 

holds for ji > Ch~^, q = 1,2; in particular, as q = 2, {I, a) = (2,0) sharp remainder 
estimate {3.7) holds. 

(Hi) Under condition {3.1) sharp remainder estimate {3.7) holds for fi > Ch~^ , g = 1,2. 
(iv) Under condition {3.3) estimate TZ < Cfi^h'^ holds with S^^{x,t) = (r = 0) and 
with arbitrarily large s. 
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4 Remarks. Generalizations 

One can generalize the results stated above. 

Remark 4.1. One can get rid of condition (0.7) by method of rescaling; then all the results 
remain the same. 

Remark 4.2. Instead of operator (0.2) one can consider operator (0.1) in its full generality: 

(i) as g > 3 no modification in conditions is needed. 

(ii) As g = 0, 1, 2 non- degeneracy conditions (0.8) and (3.1) should be modified; as r = 1 
one needs to replace |VV^| by |V(\^//)| with / = /i (r = 0); as r > 2 the modification is 
more profound because we are essentially in the matrix situation. For example if /i, . . . , 
have constant multiplicities our condition looks like 

(4.1) + ^ f,T,\ + \\/{V~\Yl f^^^)) I ^ ^0 Vn > 0, . . . , r, > 

i 3 

as (r = 0) and we need to assume that (0.7) holds. 

(iii) Smoothness conditions to g^^, Fjk should be at least (/, a) and also at least (2, 1) 
(required for reduction arguments) but for they could be even stronger to get a proper 
mollification error. 

(iv) Constant multiplicity of fj which was taken for granted is no more guaranteed and 
the results of subsection 2.3 require it. 

(v) Further, condition (2.12) should be fulfilled for n{x, h) integrated with respect to x 
and it is fulfilled automatically with ^{h) = h!^ provided {V(/2//i), . . . , V(/r//i)} has at 
least rank n — 1. 

(vi) As g = 1 and non- degeneracy condition is not fulfilled, unremovable 0(/i~^)) terms 
in the canonical form could lead to a some correction term in order to save the estimate. 
These terms can appear due to variable g^^, Fjk even if fj = const and they are unremovable 
due to the third-order resonances. 

(vii) In the case of g = and the ultrastrong magnetic field one should require that for 
each a G Z"'"'' such that | J^ji^^j ~(^j)fj\ — ^ ^ach j either dj = aj or fj has a constant 
multiphcity. 

Remark 4.3. (i) The results of this article are proven in three papers [Ivr3], [Ivr4], [Ivr5] 
where the first one is dealing with 2,3-dimensional cases and the second and the third are 
dealing with the higher dimensions. One can access them from 

http : //www. math. toronto . edu: / ivrii/Research/preprints .html 
as well as relevant talks to be viewed on a computer screen rather than printed. 

(ii) I my forthcoming papers I am planning to get rid off assumption "rank(Fjfc) = 
const" ; surely results of the case g = will be no more valid. 
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